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REPORT STATUS 


This report describes procedures for treating thermally induced structural 
vibrations by finite element analyses techniques. It is intended for readers 
familiar with finite element methods but unfamiliar with thermoelasticity theory 
and procedure. The report demonstrates the applicability of finite element 
methods to the thermal vibration problem and, as such, should serve as a guide 
to users of various discrete element structural analysis programs when treating 
problems of this type. 
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ANALYSIS OF 


THERMALLY INDUCED STRUCTURAL VIBRATIONS 

BY 

FINITE ELEMENT TECHNIQUES 


By 

James B. Mason 
Goddard Space Flight Center 


SUMMARY 


A brief review of the classical continuum approach to thermoelasticity 
problems is given. It is shown that an analogous technique can be used to re- 
duce the study of thermally induced vibrations by finite element techniques to 
that of a forced response analysis of the idealized structure excited by time 
dependent equivalent mechanical loads. 

Equivalent mechanical loads for beam and rectangular plate elements are 
presented. These time dependent loads are evaluated and employed to treat two 
sample problems for which exact solutions of the differential equations of mo- 
tion exist. The samples considered are those of a simply supported beam and a 
simply supported plate exposed to a uniform step heat input over one face. Re- 
sults of the finite element analyses are compared to the exact solutions for 
these two cases. 
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ANALYSIS OF 


THERMALLY INDUCED STRUCTURAL VIBRATIONS 

BY 

FINITE ELEMENT TECHNIQUES 


INTRODUCTION 

It is the purpose of this report to indicate the general procedures to be used 
in treating thermally induced vibrations of complex structures by finite element 
techniques . Consideration of this problem is taking on new importance in both 
present and planned spacecraft structures because of their lightweight thin 
walled construction, their associated slender structural components (typified by 
long antennas and booms) and the unusually demanding performance requirements 
expected of many of these vehicles in their space environment. 

The present investigation demonstrates that the thermal vibration problem 
can be treated by a forced response analysis of the idealized structure excited 
by time dependent equivalent mechanical loads. The problem is thus reduced to 
the determination of the equivalent mechanical loads associated with the particular 
finite elements used in the idealization of the real structure. 

Time dependent mechanical loads are obtained in this report for a beam 
element of uniform cross section and a thin rectangular plate bending element. 

In order to demonstrate the analytical procedures, two sample problems for 
which exact solutions exist are analyzed using the beam and plate elements. 

The finite element computations were performed on the IBM 7094 computer using 
the Martin Company SB-038 computer program and the results of this analysis 
compared to the exact solutions . 


ASSUMPTIONS AND BACKGROUND 

In the following development all of the usual assumptions of linear elasticity 
for continuous, homogeneous, and isotropic materials are employed and Hooke’s 
law is assumed valid. In addition, the conversion of mechanical energy into heat 
is neglected. With the latter assumption, the determination of the temperature 
distribution in the structure is uncoupled from the elasticity problem and may 
be obtained by the Fourier heat conduction equation. 

A brief review of the method of restraint for treating thermoelasticity 
problems by the techniques of classical elasticity is given here. This is done 
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to provide insight for the solution of these problems by the finite element 
approach. 

Static thermoelasticity problems may be considered as initial strain problems 
and, in turn, the initial strains treated as additional mechanical loads acting on 
the structure. The field equations describing the heated structure in rectangular 
Cartesian coordinates may be expressed as: 


Equilibrium Equation 


O'. . 

> I 


+ X; = 0 


(l) 


Stress-Strain 
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Strain-Displacement 
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where the conventional summation convention is adapted, (i, j , j3 = 1, 2, 3) 

B ( ) 

and ( ), A indicates - ^ — , here 

i 

x A ~ coordinate component 

u i ~ displacement component 

X i ~ body force component 

o\. ~ stress component 

c ^ j ~ strain component 

S ij ~ Kronecker delta 

a ~ coefficient of thermal expansion 

T ~ temperature rise from a given base temperature 

v ~ Poisson's ratio 

E ~ modulus of elasticity 
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subject to the boundary conditions: 


Traction Boundary Conditions 


n 

s. 
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Displacement Boundary Conditions 


u. 
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Mp 


(5) 


where 

n 

S. ~ surface traction component acting on the surface with outer 
normal n 

n i ~ component of n 

P (p) ~ prescribed displacement function component at points p of the 
boundary surface 

Employing the method of restraint, we consider a body free of surface trac- 
tions to be subdivided into infinitely small cubical elements. Initially these 
small elements will fit together to form a continuum. If the temperature of the 
body is altered in a nonuniform manner, however, the cubical elements will ex- 
pand an amount 


e. . 

i j 


aT S. . 
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proportional to their own temperature change and no longer necessarily form a 
continuum. This expansion can be entirely suppressed by applying a uniform 
pressure 


cr. . 

i J 


<xET 

~ 1 - 2v S ij 
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to each of the individual elements. Under this stress system, referred to as 
"state A," the lack of fit between elements in the heated structure is removed 
and the body assumes its initial undeformed continuous configuration. 



To this point we have determined the ' ‘ate A stress distribution, given by 
Equation (7), necessary to restrain the he. od structure in its undeformed con- 
figuration. The requirement now is to obtain equivalent meek oiical loads which, 
when applied to the heated structure, will produce this same restraining stress 
state. This technique of determining equivalent external loads is used in order 
that die well-known methods from classical elasticity for the analysis of struc- 
tures under external loading can be utilized in the analysis of thermally deformed 
structures. 

To obtain the equivalent mechanical loads, we note that the state A stress 
system ensures compatibility in that it is the stress system required to restrain 
the heated structure in its undeformed configuration. Substituting Equation (7) 
into the equilibrium Equation (1) results in the required body forces for restraint 


X 
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aE 
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while substitution of Equation (7) into Equation (4) results in the required re- 
straining surface tractions 


n 

s 
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aET 

1 - 2v n i 


(9) 


Therefore, Equations (8) and (9) express the equivalent mechanical loads 
which produce the state A stress distribution and restrain the heated structure 
in its undeformed configuration . 

Since in actuality these loads do not exist on the body, the undeformed struc- 
ture is now loaded with the negative of the calculated body forces 


X 

1 



( 10 ) 


and the negative of the calculated surface tractions 

“ aET 

S. - j I n. 

* 1 - 2v » 


(ID 


in a manner which satisfies the boundary conditions existing on the structure. 
The solution of this problem is obtained by classical elasticity methods for a 
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structure under the action of external loads given by Equations (10) and (11) and 
is referred to here as the "state B" solution. 

The final displacements are obtained from the state B analysis, since the 
state A displacements are zero, and the final stresses in the body are found by 
g umming the state A and state B stress distributions. 

It is seen that the effects of initial thermal strains have been treated by the 
introduction of equivalent external body and surface forces which were derived 
by restraining the deformations of the heated structure. The thermoelasticity 
problem can be treated by finite element methods in a manner completely analo- 
gous to the continuum approach described above. Conceptically, this is achieved 
by replacing the infini tesimal cubical elements and their associated restraining 
stress systems by finite elements and the necessary equivalent mechanical loads 
which must be applied at the element grid points to constrain the elements in 
their unheated shapes. 


DISCRETE ELEMENT FORMULATION 

The matrix-displacement formulation of the thermally induced vibration 
problem is used in the following development. The results, however, are appli- 
cable to other finite element formulations. 

hi the finite element approach to structural analysis problems, the body to 
be analyzed is idealized to a new structure composed of an assemblage of finite 
structural elements which are joined at grid points. The individual discrete ele- 
ments are defined by a number of grid point degrees of freedom which are suffi- 
cient to adequately represent the stress and displacement behavior of the element. 
One may write the force-displacement relations for the thermally strained i th 
element by the stiffness equation 

( X>i = [k] i {u} k T ix T } i ( 12 ) 


where 

{X'} i ~ forces acting at the grid points of element i 

[k] . ~ stiffness matrix of element i 

{u} k ~ grid point displacements of element i 

{x T ) i ~ thermally equivalent mechanical loads acting at the grid points 
of element i 
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The equivalent mechanical loads (x T > i are obtained as those required to remove 
the initial thermal strains in the discrete element. As such, the analogy with 
the restraining systems, given by Equations (8) and (9), is apparent. 

Assembly of the element matrices. Equation (12), to form the idealized 
structure results in the following set of equations describing the static behavior 
of the structure 


{X} = W {u} + 


(13) 


where 

(X) ~ external applied loads acting at the grid points of the assembled 
structure 

[K3 ~ stiffness matrix of the assembled structure 

(U> ~ grid point displacements of the assembled structure 

{Xj I ~ thermally equivalent mechanical loads acting at the grid points 
of the assembled structure 

In Equation (13), the equivalent mechanical loads acting at a grid point of the 
assembled structure are obtained as the vectorial sum, at the grid point, of the 
loads (X T } 1 of all elements joined at that grid point. 

hi many practical problems the temperature variations in the structure 
change rapidly with time and the effects of inertia can not be neglected. In these 
cases, the study of thermally induced vibrations by finite element techniques 
can be accomplished by modification of Equation (13) to Include inertia effects. 
Using the well-known principle of d'Alembert, the negatives of the inertia forces 

- lM) {U} (14) 


where 

[M] ~ mass matrix of the assembled structure 

(U) ~ grid point accelerations of the assembled structure 

are treated as additional applied loads acting on the structure. In addition, the 
equivalent mechanical loads become time dependent forcing functions. Thus, 
from Equations (13) and (14) for the case of no externally applied loading, we 
can write 
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IMG (U> + [£) {U} = - (X^t)} 


(15) 


Examination of Equation (15) shows that the problem of thermally induced 
vibrations in finite element structures is reduced to that of obtaining the dynamic 
response to a forced vibration in which the external forcing functions become the 
negatives of the equivalent mechanical loads. Upon the determination of {XjJ as 
a function of time t and imposition of boundary conditions, Equation (15) is solved 
for the unknowm grid point displacements r U (t)> of the idealized structure. 

The remaindei of this report has been written to demonstrate and evaluate 
die above analysis procedures. Two sample problems, for which solutions of the 
differential equations of motion exist, have been chosen for finite element analy- 
sis. The determination of equivalent mechanical loads for the finite elements 
employed in these samples, i.e., unifoiru beams and rectangular plates, pose no 
particular difficulty. In general, how'ever, the evaluation of mechanical loads for 
various finite elements is more involved since these loads are, in reality, func- 
tions of the assumed displacement or stress pattern for the element and are 
obtained by energy minimization procedures. 


THERMALLY INDUCED BEAM VIBRATIONS 

The problem to be considered in this section is that of a uniform simply 
supported rectangular beam of length L and thickness h exposed to a step heat 

input Q (constant in x) along edge y = + -y while edge y = - ^ is insulated, see 

Figure 1. For this case the time varying temperature distribution in the beam, 
provided T(y,o) = 0, is (1) 




where 


■© 


nondimensional time parameter 

thermal conductivity 
thermal diffusivity 


7 




Figure 1. Simply Supported Beam Exposed to Step Heat Input 


Elasticity Formulation* 

The governing equation for the displacement of a uniform Bernoulli -Euler 
beam subjected to heating can be written as 


3 2 / 3 2 v\ 3 2 v ^ 

(El + p A = 

3x 2 \ 2x 2 / 3t 2 3x 2 

where the forcing moment is 

^(t) = / f E a T (y,r) ydA 

A 


and 


x ~ longitudinal coordinate 
y ~ transverse coordinate 
v ~ transverse displacement 
I ~ area moment of inertia 
A ~ cross sectional area 
p ~ mass density 


(17) 


(18) 
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Substituting Equation (16) into Equation (18) and integrating yield the time 
dependent moment in terms of the basic input parameters as 


*t( t > 


48EIQa 

7T 4 

7T 4 k 

96 


<D 

I 




( 19 ) 


Defining the following dimensionless quantities 


£ =-r- 


OIt(t) = 


n'kMr 

192EIQa 


V(^,T) = 


A 

ff^kv 

192QaL 2 


( 20 ) 


B = 



% 


and using Equation (19), noting that M T does not depend on x, the governing equa- 
tion becomes 


. B 4 V B 2 V 

B 4 + = 0 

Bf 4 Bt 2 


( 21 ) 


The boundary and initial conditions for the simply supported beam are 


BV 


V(0,T) = V(l,r) = V(f,0) = -^r(f.O) = 0 


( 22 ) 


3*V(0,t) _ 3 2 V 

if* if* 


(1,t) - 
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Solution of the system given by Equations (21) and (22) yields the displacement 
function 


Vtf.T) 




n*l. 3, 5 


sin nrr£ 
n 


— sin n 2 w 2 B 2 r 

8B 2 n 2 


OD 

z: 


e -il*3r 


+ 



sin n 2 7r 2 B 2 r - cos n 2 7r 2 B 2 r 
j 4 + ^B 4 


(23) 


Results are given for V versus r for the case B = 1 and £ = % in Figure 2. 
The results were obtained by a single precision program which summed to 
n = j = 21 in Equation (23). 

Finite Element Formulation 

It has been shown, see Equation 12, that the stiffness equation 

(x} £ = [k]. {u^ + (24) 


represents the force -displacement relations for a thermally strained finite ele- 
ment. The stiffness and mass matrices for beams are well-known, and the 
problem of thermally induced vibrations in finite element beam problems reduces 
to the determination of equivalent mechanical loads. 

A typical beam element of uniform rectangular cross section is shown in 
Figure 3. Only temperature distributions T (y, t), constant along the axis of the 
element, need be considered if care is taken to model the structure such that 
stepwise variations in temperature closely approximate the axial temperature 
variation in the actual structure. Equivalent mechanical loads necessary to re- 
strain the beam element from bending when subjected to a temperature variation 
through the depth of the beam element are required. 

Since the temperature does not vary along the length of the element and 
equilibrium is considered in the axial direction only, no body forces are required 
for restraint and only the boundary stresses 
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NONDIMENSIONAL MIDPOINT DISPLACEMENT V 



NONDIMENSIONAL TIME r 

Figure 2. Thermally Induced Beam Vibrations (Theoretical and Finite Element Midpoint Displacement 

Response Predictions for Beam Forcing Function m^) 




need be applied. This boundary stress distribution can be replaced by the re- 
sultant equivalent mechanical moments 

-h/2 

M b T ( r ) = I baET(y,r)y dy (26) 

J- h/2 

acting at the grid points of the element, see Figure 3. 

Referring to Figure 3 and Equation (26), we see that the stiffness equation 
for a typical heated beam element can be written as: 



Having found the equivalent mechanical loads for the beam element, we can 
determine the required forcing functions for the complete model. The assembled 
finite element model of the simply supported beam of Figure 1 consisted of 11 
grid points connected by 10 beam elements of length L/10, see Figure 4. For 
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Figure 4. Finite Element Model of Simply Supported Beam 


the temperature distribution given by Equation (16), the equivalent mechanical 
moments acting on each beam element are obtained from Equation (26), assuming 
E an'* a constant, as 



48EIQa 

TT 4 k 




(28) 


Since the temperature distribution does not vary along the length in this sample 
problem, the vectorial sum of the thermal moments at the internal grid points 
(2 through 10) of the beam are zero, and the problem reduces, according to 
Equation (15), to that shown in Figure 4. 

Solution of the finite element model subjected to the time varying end 
moments given by Equation (28) and shown in Figure 4 was obtained on the 
IBM 7094 computer using the Martin Company SB-038 force method program. 
Lumped mass and modal displacement techniques were utilized for this analysis. 
Results are given for V versus r for the case B = 1 and £ = Vi in Figure 2. A 
plot of nondimensional forcing function for this case 


rr 4 k 

m b T ( r ) 192EIQa M b T 


(29) 


is also shown in this figure. 


THERMALLY INDUCED PLATE VIBRATIONS 

The problem considered in this section is that of a simply supported 
square plate subjected to a step heat input Q (constant in x and y) on face 
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z = +~2 while face z = --j is insulated, see Figure 5. The temperature distri- 
bution is identical tc that of Equation (16) with y replaced by z . 



Figure 5. Simply Supported Square Plate Exposed to Step Heat Input 


Elasticity Formulation 

The governing differential equation for the square plate is: 


DV 4 w + ph 


B 2 w 
Bt 2 



where the time dependent forcing moment is 

-h/2 

Mj(t) - | aET (z,t) zdz 

J- h/2 

and 


(30) 


( 31 ) 
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x,y ~ middle surface coordinates 
z ~ transverse coordinate 
h ~ thickness of plate 
w ~ transverse displacement 


D 


Eh 3 

12(1 - v 2 ) 


T(z,r) 



_2_ y -1 (-iy e~i 2ff2r 

772 Z-. j 2 

i = i 


cos 


j 77 
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Solution of Equation (30) for the simply supported plate yields 



r 




8 B; 





+ 


j 2 sin [B 2 77 2 (m 2 + n 2 ) r] ~| 
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where 


B 


1 



and 


a ~ square plate dimension 
Results for the following hypothetical problem 


E 


v 


P 


k 


a 


1.44 

1/3 

5.34 

1.18 

12.0 

6.90 


x 10* 


lb 

ft 3 


lb-sec 2 

ft 4 

BTU 

ft - sec- °F 


X 


10 ' 6 


ft 

ft- °F 


BTU- ft 

lb- °F- sec 2 


a - 8.0 ft 


h = 0.1 ft 


Q 


1440 


BTU 
ft 2 - sec 


(33) 
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were obtained from a single precision computer program which summed to 
m = n = j = 25 in Equation (32). Results of this calculation for midpoint dis- 
placement w versus time t are shown in Figure 6. It should be mentioned here 
that the numerical values were chosen to take advantage of an existing finite 
element plate model, see discussion of this below, and are not intended to be 
realistic. 

Finite Element Formulation 

An existing quarter symmetry model of an 8.0 ft x 8.0 ft x 0.1 ft simply 
supported plate was used in this study. The model consisted of a four by four 
network of 1.0 ft x 1.0 ft x 0.1 ft plate elements, see Figure 7. I'nlike the beam, 
however, various stiffness matrices have been suggested to describe the be- 
havior of plate elements. Hrennikoff plate elements, ( 3) obtained as an equivalent 
network of beams connected at the four corner points of the square element, 
were used. Each grid point of the model has three -degree s -of-f reedom, one 
translation (w) and two rotations (£, p). The equivalent mechanical loads acting 
at the grid points of the element were obtained by prc rating the distributed 
thermal restraining loads of a continuous square plate to the grid points of -he 
finite element. 

A rectangular plate element with an assumed temperature distribution 
T (z, t) constant over the planform of the plate is considered. As with the beam, 
this element will be satisfactory if care is taken to model the structure so that 
stepwise variations in temperature closely approximate the actual planform 
temperature. Since the temperature does not vary over the middle surface of 
the plate, body forces are not required for constraint. We may regard the plate 
as free to expand in the direction normal to the middle surface and constrain 
the x, y expansion by providing the boundary stress system 


a 


XX 


yy 


aET 

l-i/’ 


a 


*y 


0 


( 34 ) 


Thus, the plate bending may be constrained by replacing the boundary stress 
system of Equation (34) by the resultant thermal moments 

f h/2 aE 

M PT (t) = I j^7 t C 2 ’ 1 ) zdz (35) 

j - h/2 

distributed along the edges of the element. 
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Figure 6. Thermally Induced Plate Vibrations (Theoretical And Finite Element Midpoint 

Displacement Response Predictions For Plate Forcing Function M ) 

P T 


8.0 



Figure 7. Quarter Symmetry Model of Simply Supported Square Plate 


Prorating the distributed moments. Equation (35), to the grid points of the 
finite plate element results in a system of equivalent mechanical moments re- 
quired for restraint. For the square plate element of the sample problem these 
thermal moments become 


“ox**) 



(t) s 


2 


(36) 


where s is the planform dimension of the square plate element, see Figure 8. 


The stiffness equation for a typical square plate element, referring to 
Figure 8, can be written as 
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RIGHT HAND RULE 


Figure 8. Grid Point Forces, Displacements and Thermal Loads for Finite Plate Element 


20 



The equivalent mechanical loads for the 1 ft x 1 ft x 0.1 ft plate element 
with temperature distribution Equation (16) are found using Equations (35), and 
(36) to be 
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*p T (t) 


2 a QEh 3 7 t 4 

(1 - V ) 7T 4 k J** 



(38) 


when E and a are taken as constants. Since the temperature distribution does 
not vary over the planform of the plate in the sample problem, the vectorial sum 
of the thermal moments at the internal grid points of the assembled plate model 
are zero and the problem, according to Equation (15), reduces to the quarter 
symmetry model loaded as shown in Figure 9. 


SIMPLY SUPPORTED 



X 


Figure 9. Quarter Symmetry Hrennikoff Model of Simply Supported Square Plate 
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Solution of the finite element model shown in Figure 9 was obtained by the 
Martin Company’s SB- 038 program. Lumped mass and modal acceleration tech- 
niques were utilized and the results are given for midpoint displacement w 
versus time t in Figure 6. A plot of forcing functions M pT versus time is also 
shown. 


CONCLUSIONS 

It has been shown that structural vibrations induced by uncoupled transient 
temperature distributions may be treated by a forced response analysis of the 
idealized structure excited by time dependent mechanical loads applied at the 
grid points of the model. Time dependent mechanical loads for a uniform rec- 
tangular beam element and a thin rectangular plate bending element have been 
presented. 

Two thermal vibration problems have been solved by the finite element 
method: 

1. Simply supported beam subjected to a step heat input over one face: 

The finite element results compare favorably with those of the elasticity 
solution. Differences between the finite element results and the elas- 
ticity solution can be accounted for by the coarseness of partition in the 
beam model and numerical inaccuracies inherent in obtaining both 
solutions. 

2. Simply supported plate subjected to a step heat input over one face: 

The results obtained from the finite element analysis compare well with 
those from the elasticity solution. Differences can be attributed to the 
coarseness of partition in the plate model, inadequacies inherent in the 
Hrennikoff plate bending element and inaccuracies in numerical evalua- 
tion of both solutions. 

At the option of the analyst, improvements in the finite element results for both 
of the sample problems could be obtained by use of a finer element partition. 

The results of this study demonstrate that the finite element method offers 
a powerful tool for use in the analysis of thermally induced vibrations. Future 
work in this area will concentrate on the automatic calculation of equivalent 
mechanical loads for the various finite elements incorporated in the NASTRAN 
computer program. 
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